Direct URCA process in neutron stars with strong magnetic fields. 



L. B. Leinson 

Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation 
RAS, 142092 Troitsk, Moscow Region, Russia 
e-mail : leinson@charley.izmiran.rssi.ru 
Phone : (7 095) 3340912 
Fax: (7 095) 3340124 



A. Perez 

Departamento de Fisica Teorica, Universidad de Valencia 
4-6100 Burjassot (Valencia) Spain 
- - - e-mail : Armando.Perez@uv.es 

2^ ! Phone : (34 96) 3864593 

Fax: (34 96) 3642345 



0\ 



We calculate the emissivity for the direct URCA process in strongly magnetized, degenerate 
matter in neutron stars, under /3-equilibrium. We show that, if the magnetic field is large enough for 
' protons and electrons to be confined to the ground Landau levels, the field-free threshold condition 

on proton concentration no longer holds, and direct URCA reactions are open for an arbitrary 
proton concentration. Direct URCA process leads to an early phase of fast neutron star cooling. 
' This circumstance allows us to constrain the initial magnetic field inside observed pulsars. 

(N 
(N 

> • PACS : 26.60.+c,97.60.Jd 

, Keywords : Urea, Strong magnetic fields. Neutron stars 



O 



I. INTRODUCTION 



a^ ^ 

(-H \ It is now well established that pulsars are neutron stars with a very strong magnetic field. The surface magnetic 
pi • field manifests itself by synchrotron radiation from the pulsar magnetosphere and its magnitude for young neutron 
stars has been estimated to be as large as 10^^ — 10^^ G. The internal magnetic field may not manifest at the surface 
and therefore its strength is controversial but, because of the high conductivity of the core, the value of the magnetic 
field is expected to be further increased . Considering the flux throughout the star as a simply trapped primordial 
^ flux, one obtains an internal magnetic field intensity of the order 10^^ G, even if one neglects a possible twisting of 
IL* the field, an effect that would intensify the field additionally by several orders of magnitude. The time duration of 
. !^ \ the large-scale magnetic field decay via the mechanism of ambipolar diffusion might be about a few decades or a 
■ hundred year 1^. Such a simple hypothesis meets, however, a reasonable criticism, because nascent neutron stars are 
?H ■ convective. Some works have been devoted to develop an alternative to the origin of the internal mag netic field [|, 
. 5^ I based on the convective dynamo mechanism, which gives rise to a magnetic field of the order 10^^ — 10^^ G inside 
neutron stars. However, in spite of the fact that we do not know yet any appropriate mechanism to produce more 
intense magnetic fields, the scalar virial theorem allows the field magnitude to be as large as 10^^ G. Magnetic 
field of this magnitude, and even more intense, have been recently discussed Q, [Q, [|| as possibly existing 
inside some neutron stars and other compact astrophysical objects. To constrain the magnetic field inside neutron 
stars, it would be helpful to search for any observable consequences of a superstrong magnetic field. For example, one 
can observe a rapid motion of the neutron star under the recoil action, due to anisotropic neutrino emission induced 
by an intense magnetic field , or one can analyze the effect of a magnetic field on the neutron star mass . 

Moreover, if superstrong magnetic fields exist in the interior of young neutron stars, new processes would be open, 
and standard processes would be also modified by the magnetic field, thus changing the evolution of those stars. 
Studying the influence of these modifications and confrontation with observational data can be used to put some 
bounds on the initial magnetic field. The aim of this paper is to consider one of these modifications, and to motivate 
further research on this field. 

As we will show, a crucial mechanism when considering the cooling of neutron stars with superstrong magnetic 
fields are the direct URCA reactions 

p -\- e~ 11 + V (1) 
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n ^ p + e + ly (2) 

where the symbol B on the top of the arrow means that the reactions rates have to be calculated for particle states 
in the magnetic field. 

We will consider degenerate, nonrelativistic nuclear matter and ultrarelativistic degenerate electrons under /?- 
equilibrium. Matter is assumed to be totally transparent for neutrinos and antineutrinos. We also assume that at 
zero temperature the neutron star is /3-stable, but at non-zero temperature T <C ef {i = n,p, e) reactions (|l|) and (|^) 
proceed near the Fermi energies ef of the participating particles. In the field- free case, for reactions p + e~ —> n + i' 
; n ^ p + + V to take place, momentum conservation near the Fermi surfaces requires the following inequality 
among the Fermi momenta of the proton (pp), the electron (kp) and neutron {qp)- 

pp + kp > qp (3) 

Together with the charge neutrality condition, the above inequality leads to the threshold for proton concentration 
Yp = Up/ {up + Un) > 1/9, where Up and n„ are the number densities of protons and neutrons, respectively. This 
means that, in the field-free case, direct URCA reactions are strongly suppressed by Pauli blocking in the neutron-rich 
nuclear matter (assuming matter is formed only by protons, neutrons and electrons). 

In magnetized media, transversal momenta of charged particles are defined only within an accuracy Ap± ~ \feB. 
We show in the following that strong magnetic fields V eB ^ qp break down condition Eq. (|^) so that direct URCA 
processes are allowed in the degenerate n-p-e system under /3-equilibrium for an arbitrary proton concentration. We 
examine reactions (Q) and (^) in a constant magnetic field large enough so that [| B > k\l (2e) . In this case the 
Fermi energy of degenerate electrons is smaller than the energy of the first excited Landau level \Jrn? -I- 2eB . In 
other words, at zero temperature, all electrons occupy only the ground Landau band and are polarized in the direction 
opposite to the magnetic field . The number density of electrons reads 

eB f'"' „ eB 



where kp is the Fermi momentum of degenerate electrons. Following charge neutrality Up = Ue, one obtains the 
proton Fermi momentum pp to be equal to kp. Therefore, if B > kp/ (2e), then degenerate protons are also confined 
to the ground Landau state and are totally polarized in the direction parallel to the external magnetic field. We 
choose our coordinate system in such a way that B = (0, 0, B). As will be obtained through the next sections, under 
the conditions mentioned above the direct URCA process (|l|) and gives raise to a large neutrino luminosity during 
the first 10 yr. of the strongly magnetized neutron star. Such large luminosity would lead to an anomalously cold 
neutron star, a fact that can be confronted with observational data and eventually used to put some bounds on the 
initial magnetic field. 

II. CALCULATION OF THE EMISSIVITY 

The emissivity for reactions (|l|) and (|^) can be computed using standard charged-current /3-decay theory. The 
matrix element for the V-A interaction reads : 

M<x = ^ y dVW„, (r) 7^ (1 - 5^75) (r) (r) 7^ (1 - 75) vl-^ (r) (5) 

where a — + {a = —) corresponds to neutrons polarized along (in the opposite direction to) the magnetic field, (r) 
and '^e represent the neutrino and electron fields , respectively, while '^p and stand for the proton and neutron; 
G = Gpcosdc with 9c being the Cabibbo angle, and gA — 1.261 is the axial-vector coupling constant. By using 
the asymmetric (Landau) gauge A ={0, Bx,0) , the four-dimensional polarized wave functions can be expressed in 
terms of stationary states in the normalization volume V = L1L2L3. The electron wave function, corresponding to 
the ground Landau band with energy Eq = m? -\- /e| reads 



^Here e is the absolute value of electron charge. For electrons having a Fermi momentum of 50 Mev one obtains B > 2 x 10^ 
G. We set h = c = ks = 1- 
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*e (r) = —=Jj== exp{-ieot) exp [i {k^z + k2y)] foiOue (6) 

The function </3o(0 is the eigenfunction of the one-dhTicnsional harmonic oscihator, normahzed with respect to ^ = 
VTB {x + k2/ (eB)) 

D\ 1/4 



— ] expi-eBe/2) (7) 

and 
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(8) 



with fc2,^3 the electron momenta in the y and z directions, respectively. In order to take into account strong 
interactions in nuclear matter, we consider nonrelativistic protons, with an effective mass M*, moving in the self- 
consistent uniform nuclear potential Up. Then the proton wave function reads 



*P (i") = J-r- cxp{-iEot) exp [i {p3Z + p2y)] (fo (C) ( 



(9) 



where p2,P3 are proton momenta in the y and z directions, respectively, and — \J eB {x — P2/ {eB)). Polarized protons 
are given by the non-relativistic spin- up bispinor W+ normalized by the condition = 1. By introducing the 

anomalous magnetic moment of the proton, its energy reads [l2| Eq ~ M.p + £p + C/p where £p = v\l 

{2M^ and 

Mp^M;-^^{gp-\)B (10) 

with the proton's Lande factor gp = 2.79. In a similar way, we consider non-relativistic neutrons, of effective mass M*, 
moving in the self-consistent uniform nuclear potential J7„. One can describe neutron polarization by non-relativistic 
bispinors Wa {p = ±). They are normalized by the condition W\Wa = 1. Then the neutron wave function reads 



*„<T (r) = ^ exp(-z£;g„i) exp \i {q^z + q2y + qix)] ( ) (11) 
\J L2i^-iL,-i \ ^ 



Here the energy of a neutron with spin polarization a is Eq„ — M* + Sqa + Un with 

^q-y = - f^f^nB (12) 

which incorporates the neutron interaction with the external magnetic field due to its anomalous magnetic moment 
fJ-n = g-nG/ {2Mp) with gn — —1.91 . In accordance to Eq.(p^), neutrons with different polarizations u = ± have 
different momenta qF± at the Fermi surface eq± = e^. Taking into account 

= (13) 
one finds q^p_ = + AM*iinB. Explicit evaluation of the matrix elements Ma- in (^ yields 



The neutrino momentum is k = {ni, K2, K3) and its energy is lo = \k\. We have neglected the neutrino momentum 
K ^ Lo ^ T, when compared to neutron momenta qF±- The neutrino emissivity reads : 
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2tj2en 



\M+f S {E,+ +LO-Eo-ea){l- fn+) fefp 



(15) 



Here (1 — fna) fefp a-re statistical factors corresponding to Fermi-Dirac distributions of particles with chemical po- 
tentials V'n — M* + J7„ + V'e — and i/'p — Mp + Sp + Up. Fermi energy of neutrons is defined in (p^. For 
electrons and protons they are given by : 

^ 2M„ 



With these definitions one has 



fn± 
fe 
fp 
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(16) 

(17) 
(18) 

exp [{sp - e^p) IT] + 1 ^^^^ 

Since the integrand in Eq. ( [l5| ) depends on p2 and k2 only through the combination P2 + k2, we can introduce a new 
variable fc2 + P2 — >■ P2 and perform integration over k2 ■ We can also neglect the small neutrino momentum in the 
momentum conservation (5-functions. Thus, we obtain 



exp [{eg± ~ eP) /T] + 1 
1 

exp[(£o"ef)/T] + l 



G^eB 



dk3dp2dp3d^qd^K S {q^ - - k^) S (52 - P2) 



(eo + h) 



£0 



■ exp 



2eB 



i (1 + gl) {oj + K3) S (m* - Mp + Un-Up + eg+ + uj ~ Sp - eo) {1 - /„+) fefp 



+2g\ {uj - Kg) 5 (^M* - Mp + C/„ - Up + Eq- + uj ~ Sp ~ eoj (1 - /„_) fefp 
The integrals can be performed in a straightforward way by using the following identity 



(20) 



/e (eo) fp (sp) (1 - fn (e„±)) S (^M* - Mp + Un-Up + eg±+i 



d£3 

xsim: 



de2 

-00 



deife (ei) /p (£2) (1 - /„ (eg)) 

- + ea + - £2 - ei ) J (eg - £,i±) 5 (£2 - £p) ^ (ei 



eo) 



(21) 



Then, due to Pauli blocking, the main contribution to the integrals comes from £1 « £g , £2 ~ £p and £3 « £„ 
By this reason, one can substitute 5 (e'^ — £„±) 5 (£p — £p) 5 (^ef — £o)instead of S (£3 — £„±) S (£2 — £p) 6 (£1 — £0) 
Considering electrons as ultrarelativistic, we obtain that (£0 + k^) is not small only when for ^3 > 0. 
Under the /3-equilibrium assumption, one has the following relationship : 



/e (£0) fp (£p) (1 - /„ (£„±)) - (1 - /e (eo)) (1 - fp (ep)) /„ (e„±) 



(22) 



Therefore, the neutron decay process (^) gives the same energy loss rate as process (|l|), although giving final antineu- 
trinos. We finally arrive to the following total neutrino plus antineutrino emissivity: 



M* Mp B 
eo — - — - — X 

PF Mp Bo 



[(1 + gl) e (4+ - 4pl) exp 
+2g\& {ql_ - Apl) exp 



9^4 



2eB 
2eB 



4(l+5i)e((?|,+)exp - 



<1f+ 
2eB 



+ 8g\e {ql^_)exp 



2eB r 



(23) 
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Where the charge neutrality condition kp = Pf was used and ep = 2^^G^M^T^ = 1.04 x 10^'^ Tg^ (ergs/cm^s), 
with Tg the temperature in units of 10^ K and Bq = M^/e ~ 1.5 x lO^^G. Since both electrons and protons are at 
the ground Landau level, the third-component momentum conservation reads either qs = kp + PF-, when proton and 
electron go in the same direction, along the magnetic field, or = kp —pp when the initial particles move in opposite 
directions . Since kp = pp and q^^ > gf, four possibilities are open : qp_^_ > {2pp) , qp_ > {2pp) , q^^ > and 
ql_ > 0. This is the origin of the four terms in the latter equation, which are multiplied by their corresponding step 
functions, where G (x) gives +1 if the argument exceeds 0, and is otherwise. As we will show later, opening or 
closing of these channels will give rise to jumps on the emissivity. 



III. RESULTS 



In order to estimate numerically the emissivity of the direct URCA process in a strongly magnetized neutron star, 
we consider nuclear matter and electrons in /3-eQuilibrium within a nonrelativistic Hartree approach in the linear 
a-uj-p model, as it was made (relativistically) in ||8|], but with taking into account the anomalous magnetic moments 
of proton and neutron in the energies for these particles. Then the interaction energies for protons and neutrons are 
given by 

U^=f^)\, + U^Y p, (24) 
\m^J 4 \mpj 

U„=f9^)\,^U9p_\\^ (25) 
\m^J 4 \mpj 

where ni, — Up + n„ is the total number density of baryons and — Up ~ Un- Also, 

M2 B 



n 



27r2 Bo 

67r2 



irPF (26) 



^{qU + i'p-) (27) 



The parameters for the coupling constants and mesons masses are taken from [|l3| to be (M/m^)'^ = 273.87, and 
5pi 



-2 (M/nip) = 97. In this model, mass M = 939MeV is the same for proton and neutron. 



To make our estimates for the emissivity, we solved for variables qp+ and pp the /3-equilibrium condition 

V'p + = V'« (28) 

with the help of equations ( p^j27| ) and qp_ = qp^ + 4:M*finB, for a given baryon density rib. In order to be fully 
consistent with incorporating magnetic moments of proton and neutron, one would also like to find the nucleon effective 
mass M* from this modified model. However, this would add some (perhaps unessential) calculation problems. Our 
purpose here is to point out some qualitative new results and motivate more realistic computations. For this reason , 
we have chosen for our estimates the values of the effective mass that are given by . More precisely, we have taken 
M* = 0.6M for Ub equal to nuclear saturation density, which is no — 0.15/to~^ in this model, and M* = Q.3M for 
Ub = 2no. 

In Fig. 1 we have plotted the emissivity of the direct URCA process, calculated from Eq. (^Sj) as a function of the 
magnetic field intensity (in units of B/Bq) for two values of the baryon density : Ub = no and nb = 2no. We have 
chosen, as a representative value for the temperature, T = lO^K. As one can see from this figure, if the magnetic 
field is larger than i? ~ 7 x lO^^G , the I'V emissivity from the neutron star is many orders of magnitude larger than 
that in the standard cooling scenario. Another feature which is readily observed is that, as magnetic field increases, 
jumps on the emissivity can appear. The first jump in the curve = no is due to qp_ — 4p|, becoming negative 
when B/Bq > 0.014. The second, much smaller one, at B/Bo ^ 0.026 , corresponds to qp_ becoming zero. For the 
curve nb = 2no we observe a different situation : qp_ — Ap^p is always negative, but an abrupt increase appears due 
to the fact that — 4p|, is first negative and, as magnetic field increases, it becomes positive at B/Bo ~ 0.03. 

Due to the transitions discussed above, it becomes difficult to find an overall analytic formula to fit the emissivity 
dependence upon the magnetic field intensity, although it is possible to find a fitting formula which describes the 
gross tendency of e (without the details described above) for a particular value of the baryon density. We have found, 
when nb = no that the emissivity can be approximated by : 
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e ~ 2.40 X lO^^fe (2.855 - 4.2556) (^23.05e-° "^^/'' + 5.7636-° "^/'') T^erg cm-^fi-^ (29) 
where b = B/Bq . 

We now discuss the imphcations of the large emissivity produced by the URCA reactions in the presence of a strong 
magnetic field. Of course, in order to make detailed calculations, one needs to incorporate this emissivity into an 
elaborated neutron star cooling code. Aside from this, it is clear that other processes will also be modified in presence 
of the strong magnetic field (and maybe new processes will be possible). However, such considerations are out of the 
scope of this paper, as we discussed in the introduction. Here, we will only present some simple estimates. 

As it is evident from Fig. 1, the emissivity depends very strongly on the magnetic field intensity. One expects 
that an initial strong B will decay within a time Ib ■ Hence, even if the luminosity due to the URCA process is very 
high at the beginning, it will become very small for t >> ts- According to the characteristic decaying time for a 
superstrong magnetic field is ts ^ 100 yr. In order to incorporate this effect, we have taken, as the time variation of 
B : 

B{t) = B,„,t„/e-*/*« (30) 

Let us also assume that the thermal energy Uth of matter can be approximated by that of a degenerate neutron 
gas, then per unit volume one has : 

utn^^^ \at' (31) 

Since a neutron interacts with the external magnetic field only due to its anomalous magnetic moment, the thermal 
energy of a degenerate neutron gas will only slightly change under the action of the magnetic field. Neglecting this 

difference we assume A = (f )^^'^ n^^'^M*. 



One has to solve the equation 



which can be rewritten as 



dt 



dT 



dt [duth/dT) 



(32) 



(33) 



For our estimate, we will consider a fixed value of the density, with rib = Uq . Then uth ^ 6.0 x lO^^Tg erg cm~^ 
. We have also chosen an initial value of the magnetic field Binuiai = 7 x lO^^G and the temperature T(0) = lO^K. 
For the emissivity, we use the approximate formula Eq. (pof). Equation (|3^) can then be solved numerically. The 
result is plotted in Fig. 2, where the vertical axes represents Tg and the abscissa gives the time (in years). As seen 
from this figure, the temperature drops very fast at the beginning due to the URCA emission. However, as B goes 
to zero, the emission will decrease and the temperature stabilizes to a value which, in this case, is of the order of 
lO^K . In Fig. 3 we have plotted the URCA luminosity that arises from the above simplified model, with an assumed 
neutron star radius R — lO^cm. Here, one observes the large luminosity at early times and its rapid decrease due 
to the temperature drop. There is an ulterior secular decrease of the luminosity as a consequence of the magnetic 
field decay, so that LuncA will become negligeable for t » ts ■ Standard emission processes will then dominate the 
neutron star cooling. 

IV. CONCLUSIONS 

In this work we have considered the direct URCA process in degenerate nuclear matter under /3-equilibrium in a 
strong magnetic field. These reactions are strongly suppressed by Pauli blocking in the field-free case, but we showed 
that, in a magnetic field B > kp/ (2e), direct URCA reactions are open for an arbitrary proton concentration, and 
will lead to a phase of fast neutron star cooling. A simple estimate shows that, after a rapid cooling epoch (around 
10 yr. or so), the neutron star temperature becomes one order of magnitude lower than in the field- free case, a 
result which is in discrepancy with observations of pulsar data |^. One then can conclude that neutron stars with 



^For a review on pulsar data, see for example jl^ . 
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characteristic cooling times of the order 10^ — 10^ years, as in the standard scenario, could not have an initial magnetic 
field comparable or larger than the critical value Be > kp/ (2e), at least in the whole volume of the central core. 
On the other hand, if such stars exist, one has to search for anomalously cold neutron stars. In order to investigate 
the implications of this process, we have considered a very simple model. To simplify calculations we also restricted 
ourselves to the above values of the magnetic field. This approach makes possible to treat all protons and electrons 
as occupying only the ground Landau bands at zero temperature. However, as one can see from Fig. 1, the direct 
URCA process can substantially contribute to neutron star cooling, even if the magnetic field is smaller than the value 
considered above. For the URCA processes to be effective one needs the electron momentum uncertainty Afc_L ~ V eB 
to be of the order of the transversal momentum itself k±_ '--^ V2neB. This condition is still valid when electrons occupy 
a few Landau bands with n ~ 1. We also have considered matter densities not much larger than saturation densities. 
For these densities, nucleons can be treated as non-relativistic particles. More realistic calculation imply additional 
complications. If one intends to go to densities larger than the saturation density, then relativistic effects for nucleons 
should be taken into account and the many-level approximation should be used in order to calculate the energy loss 
via the direct URCA process. The standard processes should be also modified with taking into account a superstrong 
magnetic field. By making such a detailed simulation one could put more stringent bounds on the initial intensity of 
the magnetic field. 
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FIG. 1. Neutrino emissivity of the direct URCA process as a function of the magnetic field intensity B {B in units of the 
proton critical field Bq), for two values of the baryon density and for a temperature T = 10® K. 



t (yr.) 



FIG. 2. Temperature evolution, in units of 10 K for the simple model discussed on the text. 
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FIG. 3. Luminosity of the URCA process (in erg/s) as a function of time, for the same model, in logarithmic scale. 



9 



